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We investigate the production of direct photons in central Au-Au collisions at the relativistic 
Heavy-Ion Collider (RHIC) at 200 GeV per nucleon, considering all possible sources. We treat 
thermal photons emitted from a quark-gluon plasma and from a hadron gas, based on a realistic 
thermodynamic expansion. Hard photons from elementary nucleon-nucleon scatterings are included: 
primordial elementary scatterings are certainly dominant at large transverse momenta, but also 
secondary photons from jet fragmentation and jet-photon conversion cannot be ignored. In both 
cases we study the effect of energy loss, and we also consider photons emitted from bremsstrahlung 
gluons via fragmentation. 

PACS numbers: 25.75.-q,12.38.Mh 



I. INTRODUCTION 

The formation and observation of a quark-gluon 
plasma in heavy ion collisions is an important goal of 
modern nuclear physics. Even if a quasiequilibrated 
plasma is created for a brief time in the collisions, it is 
still a challenge to infer knowledge of the plasma from 
particle production. Among those proposed "probes" of 
the plasma are the directly produced real photons. Those 
photons interact only electromagnetically, and so their 
mean free paths are typically much larger than the trans- 
verse size of the hot dense matter created in the collisions. 
As a result, high-energy photons produced in the inte- 
rior of the plasma usually pass through the surrounding 
matter without interaction, carrying information directly 
from wherever they were formed to the detector. 

One expects to deduce properties of the highly excited 
matter, in particular it's space-time evolution. Direct 
photons may provide these informations in two different 
ways: via so-called "thermal photons" emitted directly 
from hot and dense matter, and via secondary photons 
from initially produced jets. There are two kinds of 
such secondary photons: they may originate from the 
"fragmentation" of jets, or the "conversion" of jets into 
photons via interactions of jets with the partons of the 
equilibrated matter. Such secondary photons compete, 
however, with the photons directly produced in an initial 
hard nucleon-nucleon scattering, in the following referred 
to as "primordial NN scattering" contribution, so a care- 
fully study of these photons is very important. In the 
following, we will shortly review these different photon 
sources. 

Primordial NN scattering. The direct photon pro- 
duction via Compton scattering and quark-antiquark an- 
nihilation can be calculated in perturbation theory, using 
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the usual parton distribution functions. In principle one 
should consider at this stage also higher order contri- 
butions, like bremsstrahlung of photons accompanying 
for example two-jet production in hards parton-parton 
scattering. However, we consider this as part of the so- 
called jet fragmentation (or bremsstrahlung) contribu- 
tion, which will be affected by the thermalize matter, 
and which we will discuss separately. 

Thermal photons. In high energy nuclear collisions, 
the density of secondary partons is so high that the 
quarks and gluons rescatter and eventually thermalize 
to form a bubble of hot quark-gluon plasma (QGP).The 
plasma expands and decreases its energy density so that 
a phase transition to hadronic gas (HG) phase appears. 
Thermal photons can be produced during the whole his- 
tory of the evolution, from the QGP phase, the mixed 
phase, and from the pure HG phase. Photons from a 
thermal source are exponentially damped so that the con- 
tribution to very high is negligible. However, its con- 
tribution to low pt is dominant. 

Jet-photon conversion. When jets pass through 
thermalized matter, they may interact. In case of the 
quark-gluon plasma, these interactions are elastic col- 
lisions between jets and deconfined partons. There is 
first of all quark-antiquark annihilation and quark-gluon 
Compton scattering, which both can produce a photon. 
These photons are affected by the plasma in two ways: 
obviously the plasma is needed to allow these interac- 
tions, but there is also a secondary effect, since the jets 
first of all lose energy in the plasma, before contributing 
to the photon production. 

Jet fragmentation or bremsstrahlung photons. 
Photon production also occurs as higher order effect in 
purely partonic initial hard scatterings: at any stage of 
the evolution of a jet (final state parton emissions), there 
is a possibility of emitting photons. Also here the pres- 
ence of a QGP will affect the results, since the jets lose 
energy during the fragmentation process. And not to for- 
get the fragmentation contribution from induced gluon 



radiation in the plasma. 

Our paper is organized as follows: in sections 2 and 
3, we discuss thermal photon production; in section 4, 
we study hard photons from primordial NN scattering 
(leading order); in section 5 we study jet-photon conver- 
sion, including the modification due to jet energy loss; in 
section 6 we compute photon production from jet frag- 
mentation, also referred to as bremsstrahlung's photons, 
again considering the effect of energy loss; in section 7 we 
finally collect our results and compare with experimental 
data. 



II. PHOTON EMISSION RATES 





Figure 1: Two-to-two particle processes contributing to the 
leading order photon emission rate. 




Thermal photon production is obtained by integrating 
the photon emission rate R (number of reactions per unit 
time per unit volume which produce a photon ) over the 
space-time history of the expanding hot and dense mat- 
ter. In this section we study the photon emission rates 
from different phases of the hot dense matter. 

The spectral photon emissivity directly reflects the dy- 
namics of real photon production reactions in thermal- 
ized matter. Commonly employed formalisms are finite- 
temperature field theory and kinetic theory. As system- 
atically studied by Kapusta et aLjH], the thermal emission 
rate of photons with energy E and momentum p from a 
small system (compared to the photon mean free path) 
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where H^'^ is the retarded photon self-energy at finite 
temperature T. This formula has been derived both per- 
turbatively and nonperturbatively. It is valid to all or- 
ders in the strong interaction. If the photon self-energy is 
approximated by carrying out a loop expansion to some 
finite order, then the formulation of Eq. |(T]) is equivalent 
to relativistic kinetic theory, where the emission rate of 
photons with energy E and momentum p from a process 
of type 1+2^3+7 reads 



^ /<n 

2—1 

X \M\- 



Spi 



{2tt)^2E, 



■){2n)H'^'^\pt+p^^-p^~p^\ 



h{Ei)f2{E2)[l±h{Ei)] 
2(27r)3 



(2) 



where the /'s are the Fermi-Dirac or Bose-Einstein dis- 
tribution functions as appropriate. Eq.|[2|) is convenient 
if the scattering amplitude, M, is evaluated in a pertur- 
bative expansion. Non-perturbative (model) calculations 
at low and intermediate energies, on the other hand, are 
more amenable to the correlator formulation, Eq.Jl]). In 
the hadronic medium, e.g., Hem can be directly related to 
vector-meson spectral functions within the vector domi- 
nance model (VDM) . Instructive investigation on photon 
emission rates from both QGP and HG phases can be 
found in [Ifnl. 



Figure 2: Bremsstrahlung and pair production contributing 
to photon emission. The bottom line in each diagram can 
represent either a quark or a gluon. 



A. Quark-Gluon Plasma 

In the thermal rate from a quark gluon plasma is 
computed using the kinetic theory formalism for the sim- 
plest two-to-two scattering diagrams such as the QCD 
Compton process qg 75 and annihilation qq ^ 57, see 
FigdJ In case of a large photon energy, and with the en- 
ergy of two initial partons being larger than the energy 
of the output photon {Ei + E2 > E :s> T), the approxi- 
mation /i(£'i)/2(£'2) exp[— (i?i + E2)/T] is employed. 
Because the light quark masses are set to zero, an in- 
frared cutoff —k^ must be placed on the four- momentum 
transfer. The infrared divergence is regulated by an in- 
finite resummation of finite-temperature Feynman dia- 
grams, following Braaten and Pisarski. This amounts to 
a careful treatment of a small part of phase space left out 
in the kinetic theory calculation by imposing the infrared 
cutoff. When the contributions from the two regions of 
phase space - below and above the cutoff - are added, 
the result is independent of the cutoff. This provides 
a parameterized thermal emission rate of photons with 
energy E and momentum p from an equilibrated QGP 
at temperature T and zero net baryon density, for large 
values oi X = E/T, given as 
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An additive "1" has been introduced in the argument of 
the logarithm to enable extrapolation to small x 

As noticed in Ref.Q , Eq.|[3l) does not yet comprise the 
full result to leading order in the strong coupling con- 
stant as- Due to collinear singularities, bremsstrahlung 
as well as pair annihilation graphs contribute at the same 
order as the resumed 2^2 processes, c.f. FiglS The 
full result, which also necessitates the incorporation of 
Landau-Pomeranchuk-Migdal (LPM) interference effects. 
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Figure 3: An interference term, involving amplitudes for 
photon emission before and after multiple scattering events, 
which contributes to the leading order emission rate. 



as shown in FiglH has been computed in Ref.[3| as 
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Figure 4: Photon production reactions and decays involv- 
ing charged pions (dashed lines) and neutral p mesons (solid 
lines). 
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We will employ the above formula, taking Nf = 3, and a 
temperature dependent running coupling constantly] 
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Effects from non-zero baryon density and from off- 
equilibrium are not included, the above rates can be con- 
sidered to be a good approximation for the midrapidity 
region of Au+Au collisions at 200 AGeV. In Figl7| we 
show photon emission rate from QGP. The daslied line 
represents the results obtained by Kapusta et a/||l|, sum- 
marized in eq.|[3]), the solid Hne refers to the calculations 
of AMY|3j|, given in eq. ^ , with additional contributions 
compared to [l|. We can see the full contribution to pho- 
ton emission from a QGP is much higher than the one 
from 2 — > 2 partonic processes. 



B. Hadronic Matter 

Photons can also be produced in a hadronic phase, 
from several elementary interactions, see FiglH The 
dominant contribution comes from the reactions tttt 
P7 and irp irj; The decay p tt+tt+j also contributes 
significantly. Interactions involving strange mesons or 
baryons can also produce photons, but these contribu- 
tions are relatively small because of the phase space sup- 
pression due to their big masses. The situation of ther- 
mal photon radiation rates from a hadronic gas is un- 
certain, due to difficulties related to the strong coupling 



Figure 5: Feynman diagram of TTp ■K'y through ai(1260) 
resonance. 



and the masses of hadrons. The study is usually carried 
out within effective Lagrangians. Constraints on the in- 
teraction vertices can, to a certain extent, be imposed 
by symmetry principles (e.g., e.m. gauge and chiral in- 
variance). CoupHng constants are estimated by adjust- 
ing to measured decay branchings in the vacuum. Thus, 
for the temperature ranges relevant to practical applica- 
tions, T=100-200 MeV, the predicted emission rates are 
inevitably beset with significant uncertainties, and there- 
fore a careful judgment of the latter becomes mandatory. 

Investigations along these lines were initiated in 
Ref.[l|, where the photon self-energy has been computed 
to 2-loop order for a mesonic system consisting of sharp 
(zero width) tt-, it- and p-mesons (plus direct uj — > 7r°7 
decays). In Ref.pl it was pointed out that np — > 7r7 
scattering via ai(1260) resonance formation (or, equiv- 
alently, ai tt^ decay), c.f. FiglH constitutes an im- 
portant contribution. This was followed up by a sys- 
tematic treatment!^ of an interacting irpai system to 2- 
loop order within the Massive Yang-Mills (MYM) frame- 
work of introducing axial-/vector mesons into a chiral 
Lagrangian, and, later, within the Hidden-Local Symme- 
try (HLS) approach^. The effect on in- medium vector 
and axial-vector meson masses is studied by Song and 
Fai Q. 

We will use the results of the MYM calculation 
where photon production from strangeness bearing 
mesons has been included as well as the axial meson ai as 
exchange particle for non-strange initial states. We show 
these rates in figlH and we list the corresponding param- 
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eterized rates, in unit of fm "^GeV ^ , with the photon 
energy (E) and the temperature (T) both in GeV: 
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factor of 3 at photon energies around E ~ 2.5 GeV. 
The reduction of the rate introduced by Rapp et al. [§| in 
2 ~ 3 GeV region of photon energies amount to a factor 
of ~ 4, confirming roughly Kapusta's findings. 

In ref. ^2|], for each hadronic vertex appearing in the 
ampHtudes, it is proposed to take hadronic form factor 
for t-channel meson X-exchange according to 
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with A=lGeV and t being the average momentum trans- 
fer t — —2Emx, where E is the photon energy and mx 
the mass of the hadron X. In our calculation, we fol- 
low this procedure: The photon emission rate is taken to 
be the sum of all terms in eqs. (8-14), with each term 
multiplied by F'^{t). The corresponding curve is shown 
in Figlll as the dashed-dotted fine. The dotted fine is 
the photon emission rate without considering hadronic 
form factors. The form factors of Ec|.l|15p makes the sup- 
pression stronger compared to P, For example, at a 
photon energy oi E = 2.5 GeV, the suppression factor is 
about 10 for 7r-exchange contributions, and even 625 for 
if-exchange processes. The suppression are even larger 
for higher photon energies. So the hadronic emission rate 
is reduced considerably after including the hadronic form 
factor as in Eq. lfTSi) . The HG rate is therefore much 
smaller than QGP one (full line). 
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Parameterisations for K* K+tt+j and K+K — + p+j 
do not appear because their rates have been found to be 
neghgible. 

Hadrons are composite objects, so they may need ver- 
tex form factors to simulate finite hadronic size effect, in 
particular at high momentum transfer. How much will 
form factor influence the results? Form factors are a very 
delicate subject, especially when electromagnetism and 
its gauge invariance are involved. For those nonstrange 
reaction channels as originally studied by Kapusta et al. 
[H, considering form factors provide a typical net sup- 
pression compared to the bare graphs by an appreciable 



III. THERMAL PHOTONS FROM THE 
EXPANDING HOT AND DENSE MATTER 

The expanding thermaHzed matter is treated by em- 
ploying three-dimensional hydrodynamics. We compute 
initial conditions at some given proper time tq, expressed 
via energy density e(To), net flavor density /^(to), and 
collective velocity v{tq), by employing the EPOS model 
[3, The hydrodynamic evolution is reaHzed using 
SPheRIO [iBl, which is a "Smoothed Particle Hydrody- 
namics" implementation, a method originally developed 
in astrophysics, and later adapted to relativistic heavy 
ion collisions . The three-dimensional hydrodynamics de- 
scribes the space-time evolution of the hot dense matter 
created in heavy ion collisions, via the 3-velocity v, the 
energy density e, the entropy density s and the baryon 
number density ns, as functions of the space-time posi- 
tion (77, T, r, (j>), with r] being the space-time rapidity, and 
r, (j) the transverse coordinates. In this paper, we con- 
sider central AuAu coUisions (10% most central events) 
at 200 AGeV. The corresponding results of a hydrody- 
namical evolutions is shown in FigUl where we plot the 
r-dependence of e at = for different values of and 
T. The initial time is tq — 0.5fm/c. The solid lines 
and the corresponding dotted lines refer to 4> = and 
(f) = 7r/2, respectively. The three horizontal dotted lines 
are the the energy densities £1 = 1.675 GeV/fm'^ and 
£2 = 0.325 GeV/fm'^ , limiting the mixed phase, and the 
freeze-out energy density £3 = 0.08 GeV/fm"^.. 
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Figure 6: Photon emission rates from a HG from the different reactions (bare graphs, no hadronic form factor). 
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Figure 7: Photon emission rate from QGP and from HG. 
The dashed line represents Kapusta's et aZpj] calculations of 
photon emission rates from a QGP at 200 MeV via 2—^2 
processes, i.e., QCD Compton process and quark-antiquark 
annihilation. The solid line refers to AMY[_3], where in ad- 
dition the LPM contribution was taken into account. The 
dotted line is a sum of photon emission rates from a HG via 
different reactions (bare graphs) shown in fig. [6l The dot- 
ted dashed line stands for photon emission rates from a HG 
including hadronic form factors. 



S2 being the entropy densities corresponding to Si and £2, 
s is the total entropy, which is simply a Hnear function 
of £ , like s = ae + b. Then /qgp and /hg are linear in e 
as well, namely 



/qgp — (1 - /hg) — 



£2 



£1 - £2 



(16) 



The photon emission rates discussed in the previous 
chapter are obtained in the local rest frame, so we should 
better write 

where quantities viewed in the local rest frame are dec- 
orated with the superscript *. Bare symbols (without 
*) refer to the laboratory frame. We need the photon 
spectrum observed in the laboratory, which is given as 

dN ^ ^dN^ [d^xTiE T) 
dyd'^PT <Pp J ' 

j d^xT{E*,T) 
= J dTdxdydriTT{E*,T), (18) 



The relation between energy density and temperature, 
T = T{e), as used in SPheRIO is shown in FigE] (and 
tabulated for later use). 

We define /qgp to be the fraction of matter in the 
QGP phase and /hg as the corresponding fraction in the 
HG phase, at each space-time point {x,y,f],T). We have 
obviously /qgp = 1, /hg = 0, if the energy density is 
bigger than £1, /qgp = 0, /hg = 1, if the energy density 
is between £3 and £2, and /qgp = 0, /hg = 0, if the 
energy density is smaher than £3. In the mixed phase 
(£2 < £< £1, we have si/qgp + S2/hg = s, with si and 



with 

r{E* , T) - /QGpr'^^P^^ (E* , T) + /Hor^G-*^ {e* , t) 

(19) 

where T'^gp^^, ^ ^^^^ pHG~,7 (E*,T) are the photon 
emission rates from QGP phase and from HG phase as 
discussed in the previous chapter. The center-of-mass 
energy E* in eq. (fT8| is related to the photon momenta 
in the observer frame (appearing on the l.h.s. of Eq. (fT8|) ) 
as 

E* ='-fE-jv(T,x,y,'q)-p, (20) 
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Figure 8: The energy density of the hot dense matter created in the most 10% central Au-Au collisions at 200 AGeV, is plotted 
as function of the transverse space variable r. Curves from up top to bottom correspond to different proper times, t=0.5, 2.5, 
... (fm/c). The initial time is tq = 0.5fm/c. The solid lines and the corresponding dotted lines refer to (j) — Q and 4> = ""/S, 
respectively. The two upper horizontal dotted define the mixed phase, the lower horizontal line the freeze-out energy density 
£3 = 0.08 GeV/fm^ 
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Figure 9: The relation between temperature and energy 
density. A first-order phase transition happens when 
e{x,y,rj,T) G (e2, ei) with space-time point {x,y,rj,T) in a 
mixed phase. 



with 
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^Jl - \v{T,x,y,ri)\'^ 



(21) 



where w(r, x, y, r]) is the flow velocity at a given space- 
time point. The flow rapidity y = ln(l -I- Wz)/(1 — Vz) 
is roughly equal to the space-time rapidity (as in the 
Bjorken model). The radial dependence of the transverse 
velocity Vr is shown in Figfim for central Au-Au collisions 
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Figure 10: Transverse velocity Vr of the hydrodynamic expan- 
sion for the hot dense matter created in the most 10% central 
Au-Au collisions at 200 AGeV. 



at 200 AGeV. The large transverse flows at large radii will 
boost thermal photons to higher pt regions. 

In FigdU we present the transverse momentum spec- 
tra of thermal photons produced in the 10% most cen- 
tral Au-Au collisions at 200 AGeV with an initial time 
of To=0.5 fm/c. The contributions from the two phases 



7 



> 
o 

-d 



1 














□ ther. w/o FF 


-4 


\ 




• ther. w/FF 






* ' ■. 




-8 

10 








-12 








10 




QGP^ Y 








HG^ y w/ FF 




-16 








10 




HG^ y w/o FF 






1 1 


1 1 1 1 1 1 





10 



P, (GeV) 



Figure 11: Thermal photons from the most 10% central Au- 
Au collisions at 200 AGeV with initial time to=0.5 fm/c. The 
contributions from the two phases are presented separately: 
QGP phase (solid line) and HG phase (dashed line). The 
upper limit of the contribution from HG phase, i.e., with- 
out consideration of hadron form factor, is presented as dot- 
ted line. The total thermal contribution is plotted in the 
two cases, with hadron form factor (full circles ) and without 
hadron form factor (empty squares). 



are presented separately: QGP phase (solid line) and HG 
phase (dashed Hne). We can see thermal contribution is 
dominated by photons from QGP phase, which is much 
bigger than the upper limit of the contribution from HG 
phase, i.e., without consideration of hadron form factors 
(dotted line). The reason is simply that the QGP phase is 
"hotter" compared to the HG phase. Therefore, hadron 
form factors makes very little difference concerning the 
total thermal contribution, as seen in FiglHl when com- 
paring the complete results with form factor (full circles) 
and without form factor (empty squares). 



IV. HARD PHOTONS FROM PRIMORDIAL 
N-N SCATTERING, INCLUDING SHADOWING 
AND ISOSPIN MIXING 

Here we present the leading order perturbative calcula- 
tion of hard photon production, namely the hard photons 
from primordial N-N scattering. Higher order contribu- 
tions are related to the production and propagation of 
jets, so we treat them later. The spectrum of hard pho- 
tons from a collision between nuclei A and nuclei B is 
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where Ga/A{^a^ M'^) and Gh/sixb, M'^) are parton distri- 
bution functions for nuclei AandB. We take MRST 2001 



LO parton distributions for protons[l^. Nuclear shad- 
owing effects are taken into account by using EKS98 scale 
dependent nuclear ratios -Ra^®(a;, ^)[l§]. The mixed 
isospin in nuclei with mass A, neutron number N and 
proton number Z is taken into account as 

AT 7 

G^Ia(x) = A\-G^lj,{x)^-G^lp{x)\Rf^^(x,A). (23) 

The elementary cross sections after color sum and spin 
average are given as 1(13] 
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and Aqcd = 200 MeV. The 4-momenta of the incoming 
particles {pa, Pb) and of the photon (pj) in the center-of- 
mass are 



Pa = (a^a^, 0, 0, Xa^), Pb = i^b^, 0, 0, -Xb^) 



and 

Pj = (pTCOsh?/,pT,pTsinh?/). 
The Mandelstam variables are then given as 

S = (Pa +Pbf' XaXbS 
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u= {pb- P-rf = ~Xby/spT exp(y). 

We set the factorization scale M and renormalization 
scale Q to be pT- Due to the (^-function in Eq. l(22|) . the 
Xfc-integration is trivial, leading to 

^ XaX± exp{-y) 
2xa - x± exp(2/) ' 

with x± — 2pt/\/s. The condition Xb < I requires Xa G 
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lin, 1) with a;inin = x± exp{y)/{2 - xi_ exp(-?/)). So 
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In FigfT2l we plot the corresponding spectrum, for cen- 
tral Au-Au collisions at 200 AGeV. Empty circles present 
PHENIX datai25:]. At high pt region, the contribution 
of primordial NN collisions is very close the experimental 
data. 
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Figure 12: The contribution of primordial NN collisions to 
direct photon production in the 10% most central Au-Au col- 
lisions at 200 AGeV. Isospin mixing and nuclear shadowing 
are considered. Empty circles present PHENIX data[25l]. 



JET PHOTON CONVERSION, INCLUDING 
JET ENERGY LOSS 



When jets propagate in the hot dense matter created 
in heavy ion coUisions, they also interact with the mat- 
ter and produce direct photons via the Compton process 
19 ~^ 11 and the quark-antiquark annihilation process 
QQ ^ 91- We denote the momenta of the jet, the thermal 
parton, and the photon by pjot, Pth and p, respectively. 
The leading order QCD Compton and annihilation cross 
sections are peaked in the forward and backward direc- 
tions. In the laboratory frame we have |pjct| ^ bth| ~ T, 
where T is the temperature of the plasma. For high en- 
ergy photons, i.e. |p| ^ T , this implies that \p\ ~ bjetl- 
That is why we call this process a conversion of a jet into 
a photon with similar momentum. The rate of photon 
production by annihilation and Compton scattering of 
jets in the medium can be approximated as [20|(2l| 



In — Ir- + C 



(34) 



with C = -1.916, TOjj^ ~ and where as — 

g^/(47r) and a are the strong and the electromagnetic 
couplings. Th subscript "g" denotes all light quark and 
antiquark species with charge Cq , and fq = fq {p, x) is the 
phase-space density of partons of flavor q. It is worth em- 
phasizing that the conversion property of the process is 
reflected in eq. l(34|) by the fact that the photon spectrum 
is directly proportional to the parton spectrum fq. Then 
at midrapidity the jet-photon conversion contribution to 
the direct photon production in the most 10% Au-Au col- 
lisions is gained via the integration over the space-time 



evolution of the hot dense matter in QGP phase: 

^jyjet+QGP^7 



dy(PpT 



(35) 



E 



dN 
(Ppd^x 



fQGp{x, y, 1], T)dxdydriTdT. 



To get fq, the phase space densities of quarks and an- 
tiquarks, we have to flx the geometry of jet formation. 
The jet production from primordial N-N scattering is as- 
sumed to happen at the same proper time t — 0. Then 
at r = 0, the phase space distribution of partons of type 
q is 



fq.iP^r^T^ 0) cx Ta(x- 



-,y)TB{x 



i^.y)5{z), (36) 



where b in the impact parameter and Ta and Tb are 
thickness functions of nuclei A and B. The (5-function 
is motivated by the strong Lorentz contraction of the 
colliding nuclei. To have a simple form for a central A-A 
collision, i.e. most 10% central Au-Au collision, we take 
the approximation 



/, (p, f , T - 0) - ^-^ 3:35:^ (37) 



E dyd'^pT 



where i?j_is the radius of nuclei Au, and with r± = 
\J + y^. The momentum distribution is calculated as 



dyd'^PT 



(38) 



J2 j dXaKj,tGa/A{Xa,M^)Gb/B{xb,M^) 

abed 

XaS- \/spT exp(y) tt dt 



which is very similar to eq. l(33|) . Here, the cross sec- 
tions of all possible partonic processes ab cd make 
totally 127 terms[l3|- Kjet='2 is used to take into ac- 
count higher order contributions to jet production in our 
calculation [221 • The phase space distribution of jet sat- 
isfies the normalization 



fqiP,r,T = 0) 



d^rd^p 
(27r)3 



N, 



(39) 



where Nq is the number of jets of type q. 

If we ignore the jet energy loss due to the interaction 
between jets and matter, then at any t > 0, the phase 
space distribution of jets is 



f{p,r,r) 



d'^rQf{p,ro,T = 0)(5(r - fo - vt) 



f{p,r-vt,0) (40) 
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where 17 is the velocity ofajet, v =p/E, andt = Tcoshirf) 
with -q being the space time rapidity. 

If we consider the modification of the jet energy and its 
momentum due to the interaction between jets and hot 
dense matter, then the phase space distribution /(p, r, r) 
of jets should be replaced by 



dVod%/(po, fo, {))6{ph ~ j-AE - p)S{f- fo - vt) 



E 



= I d^PofiPo, r- vt, 0)S{p'o - ^AE - 



(41) 



with E = \p\, and where AE is the energy loss of a jet 
propagates from the formation point (ro,0) to the jet- 
medium interaction point (r, t). In static matter, one 
has fH 



AEc = ^ec 



(42) 



with Cc = as fi^E/X^, and where the mean free path Ac 
of the jet in the medium, and fj. = (?T, are all temperature 
dependent quantities. The index "c" refers to the jet type 
(quark or gluon). Due to the space-time evolution of 
temperature as T = T(a;, y, ry, r), we have to replace L 
in Eq. l(42l) by the integration v fqcpdi' along the jet's 
trajectory in the QGP plasma. Light-flavour quarks and 
gluons are massless, with v = 1, and so we get 



AE 



£c(T(x,y,7?,T)) 
>^c{T{x,y,r],T)) 



fQGp{x,y,r],T)dt'. (43) 



The mean free path Ac is given as 



jPq+aqgPg, 



(44) 



(45) 



where the cross sections are given as [24] Ci = Q ^ , with 



1 and Co 



qq — Q, Cqg = 1 and Cgg = |, and where pq and pg are 
the thermal parton densities. 

Figim shows the jet-photon conversion contribution 
to direct photon production in central Au-Au collisions 
at 200 AGeV. We show the results with and without 
energy loss, together with PHENIX datajl^l- We can 
see that considering energy loss can indeed suppress the 
jet-photon conversion contribution, by about a factor of 
three. Obviously, higher jet energies are needed to pro- 
duce a photon with a given energy, in case of jet en- 
ergy loss in the medium. The effect from different hy- 
drodynamic initial times tq is also considered. However, 
results with ro=0.5 fm/c and 1 fm/c are indistinguish- 
able. Our results also agree well with the earlier work by 
Turbide[l3|, although they use one dimensional hydrody- 
namics. 
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Figure 13: Jet-photon conversion contribution to direct pho- 
ton production in the 10% most central Au-Au collisions at 
200 AGeV. Dotted line: jets traverse the hot dense matter 
without energy loss; solid line: jets traverse the hot dense 
matter with energy loss; Empty circles present PHENIX 
data(25l|. The effect from different hydrodynamic initial times 
To is also considered. However, results with ro=0.5 fm/c and 
1 fm/c are indistinguishable. 



VI. BREMSSTRAHLUNG PHOTONS, 
INCLUDING JET ENERGY LOSS 

We discussed earlier leading order photon production 
from elementary NN collisions. However, also higher 
order diagrams give important contributions. Here, 
we discuss "bremsstrahlung photons" from jets, also re- 
ferred to as "jet fragmentation". We can treat the 
bremsstrahlung contribution via parton fragmentation 
functions D^/^{z, Q^) being the probability for obtain- 
ing a photon from a jet c, where the photon carries a 
fraction z of the jet's momentum. The effective fragmen- 
tation functions for obtaining photons from partons can 
be calculated. The leading order result is 



zD,,q{z, Q') ^el-[l + {l- zf] ln(QVA2) (46) 



and 



zD^/giz,Q^)^0, 



(47) 



where is the fractional charge of the quark q. The 
photon fragmentation functions evolve with just as 
the usual hadronic fragmentation functions do, as a result 
of gluon bremsstrahlung and qq pair production. The 
resulting evolution equations are 



dD^,qSz,Q^) 

dt 



(48) 



27r 



/ ^[D,^q^iy,Q')Pqq{z/y) 

J z y 



+D,/g{z,Q^)Pgq{z/y% 



10 



and 



dD^/g{z,( 



dt 



(49) 



2-K 



2Nf 



-[Y,D,,,^{y,Q^)P,g{z/y) 



+D^^g{z,Q')Pgg{z/y)] 

where Pqq, Pgq,Pqg and Pqq are splitting functions ap- 
pearing in the DGLAP equations. The parameterized 
solutions from Owens [ITI are 



zD^/q.{z,Q ) 

„2 ' 



(50) 



27r' * 1 - 1.63 ln(l - z) 
+0.002(1 - z)2z-i-54] ln(QVA')> 



and 



zD^/g{z, Q') = —0.0243(1 - z)z-'>-^^ ln(Q7A^). (51) 

The bremsstrahlung contribution to direct photon pro- 
duction is then 



dyd^PT 



c=g,qi 



'^'^ dyd^p^^ ;|^./c(^c,Q^), 



(52) 

with = pt/zc being the momentum carried by the jet 
c before fragmentation, and d^p/E = z'^d^p'^/E'^. The 
jet cross section can be obtained from Eq. l|38p . 

To take into account the energy loss of the jet in the 
QGP phase, the bremsstrahlung contribution to direct 
photon production is modified. We may use a modified 
fragmentation function[2^, given as 



D^/,{z„Q^,AE,) 



(53) 



= (1 - e-^)[^i?!J/,(4, Q') + ^'lD'^^^{z'g,Q')] 

Ac 



+e--i?°/,(z„Q2), 

where z'^ = Pt/{pt - ^Ec) and z'g = {L/Xc)pt/AEc 
are the rescaled momentum fractions carried by jet c and 
the emitted gluons before fragmentation, and where Z?" 

and D°^g are the original fragmentation functions, given 



in Eqs. ((5Ql) and ijSTI) . So a parton c has the probabil- 
ity exp(— L/Ac) to fragment directly without interacting 
with the medium, and the probability 1 — exp(— L/Ac) 
to interact with the medium before fragmentation. The 
ratio L/Ac represents the number of scatterings, and a 
gluon of energy Cc is emitted each time when a parton 
scatters with the medium. 

As discussed earlier, the energy loss in the medium 
with a constant temperature is AEc = j^Ec, where 



Au+Au y Vs^^=200GeV 
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20 
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Figure 14: Bremsstrahlungs (= jet fragmentation ) contri- 
bution to direct photon production in the 10% most central 
Au-Au collisions at 200 AGeV. Dotted line: jets traverse the 
hot dense matter directly without energy loss. Solid line: jets 
traverse the hot dense matter with energy loss. Empty circles: 
PHENIX datadil. 



Ac is the mean f ree path of jet c in the medium and 
Ec = OLs\J [i^E jXc. Those quantities are temperature de- 
pendent. In our case, the temperature evolves with space 
and time. So similar to Eq. l(43l) . we replace i/Ac and 
AEc by the corresponding mean values, namely 



1 



K{T{x,y,T],T)) 



ec(T(x,y,?7,r)) 
>^c{T{x,y,ri,T)) 



fQGp{x,y,ri,T)dt, (54) 



fQGp{x,y,V,T)dt. (55) 



We take the mean energy loss per scattering XAEc/L as 
the energy carried by each emitted gluon. 

In FigfT4t we present the bremsstrahlung contribution 
to direct photon production, with and without consider- 
ing energy loss. We can see that the photon production is 
considerably suppressed due to the energy loss of the jets 
in the hot dense matter. The results from two different 
initial times to=0.5 fm/c and 1 fm/c is indistinguishable. 



VII. RESULTS AND DISCUSSION 

In the following we are going to collect and discuss our 
results. In Fig[T5]we compare the different non-thermal 
contributions to photon production in central Au-Au col- 
lisions at 200 AGeV. We show the contribution from pri- 
mordial NN coUisions (dashed line), bremsstrahlung (or 
fragmentation) photons (dotted line), and photons from 
jet-photon conversion (solid line). We have taken into 
account energy loss of jets in the hot and dense matter, 
both for bremsstrahlung photons andjet-photon conver- 
sion. We also plot PHENIX data[25|, as a reference. 
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Figure 15: Comparison of the different non-thermal con- 
tributions photon production in the 10% most central Au- 
Au collisions at 200 AGeV. Dashed line: contribution from 
primordial NN collisions. Dotted line: contribution from 
Bremsstrahlung. Solid line: jet-photon conversion contribu- 
tion. Jets energy loss in the hot dense matter is considered 
both for bremsstrahlung photons and jet-photon conversion. 
Empty circles: PHENIX dataQ. 



Jet-photon conversion contributes the same magnitude 
as bremsstrahlung, but both are small compared to the 
hard photons from primordial NN collisions. 

In FigdU the total contribution (solid line), including 
all non-thermal and thermal contributions, is compared 
with PHENIX datapHENIX data] (empty circles). The 
thermal contribution from a hydrodynamic calculation 
with an initial time to=0.5 fm/c is plotted (dotted Hne) 
as well as the results of a calculation with ro = l fm/c 
(solid line with full circles). Smaller tq gives a some- 
what broader distribution, but this effect is invisible in 
the total contribution. At very low px , thermal pro- 
duction dominates the direct photon production. How- 
ever, the thermal spectra decrease very fast with pr, 
and become neghgible for pT> 4GeV/c. The hydrody- 
namic configuration also affects the jet-photon conver- 
sion and the bremsstrahlung contribution. But this ef- 
fect is very weak: taking different initial times within the 
range [0.5,l]fm/c make very little difference to the total 
contribution (invisible in Fig. 16). The complete contri- 
bution, containing all the thermal and non-thermal con- 
tributions, follow quite well the PHENIX data. 

It is useful to compare the AuAu results with proton- 
proton, and here it is convenient to study the nuclear 
modification factor, defined to be ratio of the nuclear 
spectrum to the proton-proton one, divided by the num- 
ber iVcoii of binary coUisions. We use A^coii=880 for the 
10% most central Au-Au collisions at 200 AGeV, treated 
in this paper. We compute the proton-proton differ- 
ential cross section by using (jf^^i^stic =40.83mb, and 
Kjet='2, and the unmodified photon fragmentation func- 
tions given in eq.|| 



Figure 16: The total contribution (solid line) is compared 
with PHENIX datafii] (empty circles). Thermal contri- 
butions from a hydrodynamic calculation with initial time 
ro=0.5 fm/c is plotted as dotted line, while ro = l fm/c is pre- 
sented as solid line with full circles. Different values of tq are 
indistinguishable concerning the total contribution. 
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Figure 17: Different non-thermal contributions and thermal 
contribution to nuclear modification factor Raa is compared 
with PHENIX datajii] (full squares). Dashed line: the pri- 
mordial NN contribution; Dotted line: bremsstrahlung; Solid 
line: jet-gamma conversion; Dotted dashed line: thermal con- 
tribution. 



In FigfTTl we plot the different non-thermal con- 
tributions, namely the primordial NN contribution 
(dashed line), bremsstrahlung (dotted line), and jet- 
gamma conversion (solid line). The contributions from 
bremsstrahlung and jet-photon conversion are of the 
same magnitude, but one order smaller than that the 
primordial NN curve. We also show the rapidly falling 
thermal contribution (dotted dashed line). Above 5 GeV, 
only photons from primordial NN scattering contribute 
significantly. This is why we look more closely to this 
latter contribution. 
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Figure 18: A close look at primordial NN scattering con- 
tribute to the nuclear modification factor Raa- Solid line: 
the complete primordial NN contributions, including isospin 
mixing and shadowing; Dotted line: omitting isospin mixing 
(but considering shadowing); Dashed line: omitting shadow- 
ing (but considering isospin mixing). Full squares: the same 
data as in Figll7l 



In Figim we show again the complete primordial NN 
contributions, including isospin mixing and shadowing. 
But we also show the results one would obtain by omit- 
ting isospin mixing (but considering shadowing), and 
by omitting shadowing (but considering isospin mixing) . 
Obviously, both effect are crucial. In particular the 
isospin mixing is responsible for getting a nuclear modi- 
fication factor of less than unity. 



Finally we show in Fig[T9]the complete calculation, for 
two different options of the initial time for the hydrody- 
namical evolution. The two curves are almost identical. 
As mentioned earlier, the total contribution to the nu- 
clear modification factor is less than unity, due to isospin 
correction. However, the experimental data drop sys- 
tematically below the theoretical curve, and it seems dif- 
ficult with the processes discussed in this paper, to get 
the theoretical curve further down. So maybe some "new 
physics"? Before answering this question, one should not 
forget the large uncertainties of the experimental pp ref- 
erence, which is a fit function passing between strongly 
fiuctuating data points at large pt. 
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Figure 19: The complete calculation of the nuclear modifi- 
cation factor is compared to PHENIX data(25l|. Dotted line: 
ro=0.5 fm/c; solid line: ro = l fm/c. Full squares: the same 
data as in Fig ll7l 
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